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In this paper we investigate from first principles the effect of the electron-phonon interaction 
in two paradigmatic nanostructures: irans-polyacetylene and polyethylene. We found that the 
strong electron-phonon interaction leads to the appearance of complex structures in the frequency 
dependent electronic self-energy. Those structures rule out any quasi-particle picture, and make the 
adiabatic and static approximations commonly used in the well-established Heine Allen Cardona 
(HAC) approach inadequate. We propose, instead, a fully ab-initio dynamical formulation of the 
problem within the many-body perturbation theory framework. The present dynamical theory 
reveals that the structures appearing in the electronic self-energy are connected to the existence 
of packets of correlated electron/phonon states. These states appear in the spectral functions even 
at T = QK, revealing the key role played by the zero point motion effect. We give a physical 
interpretation of these states by disclosing their internal composition by mapping the many body 
problem to the solution of an eigenvalue problem. 
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I. INTRODUCTION 

The electron-phonon (EP) coupling is well known to 
play a key role in several physical phenomena. For exam- 
ple it affects the renormalization of the electronic bandsi, 
the carriers mobility in organic devices^ or the position 
and intensity of Raman peaks^. The EP coupling is also 
the driving force that causes cxcitons dissociation at the 
donor/acceptor interface in organic photovoltaic^ and the 
transition to a superconducting phase in solids^. 

Despite the development of more powerful and efficient 
computational resources the calculation of the effects in- 
duced by the EP coupling in realistic materials remains 
a challenging task. In addition to the numerical diffi- 
culties, it has been assumed, for a long time, that this 
interaction can yield only minor corrections (of the or- 
der of meV) to the electronic levels. As a consequence 
the majority of the ab-initio simulations of the electronic 
and optical properties of a wide class of materials are 
generally performed by keeping the atoms frozen in their 
crystallographic positions. It is actually well-known that 
phonons are atomic vibrations and, as a such, can be eas- 
ily populated by increasing the temperature. This naive 
observation is de-facto used to associate the effect of the 
EP coupling to a temperature effect that vanishes as the 
temperature goes to zero. However this is not correct 
as the atoms posses an intrinsic spatial indctermination 
due to their quantum nature, that is independent on the 
temperature. These quantistic oscillations are taken into 
account by the EP coupling when T — >■ in the shape of 
a zero-point -motion effect. 

Many years ago^ Heine, Allen and Cardona (HAC) 



pointed out the EP coupling can induce corrections of 
the electronic levels as large as those induced by the 
electronic correlation. As a consequence the generally 
accepted statement that the EP coupling always yields 
minor corrections was doomed to fail. Nowadays, the 
advent of more refined numerical techniques, has made 
possible to ground the HAC approach in a fully ab- 
initio framework. This has been used to compute the 
gap renormalization in carbon-nanotubes^, the finite 
temperature optical properties of semiconductors and 
insulators^, and to confirm a large zero-point renormal- 
ization (615 meV) of the band-gap of bulk diamond^, 
previously calculated by Zollner using semi-empirical 
mcthodsi^. These works are calling into question decades 
of results, by instilling the doubt that a solely electronic 
theory may be inadequate. 

In this work we show that in nano-structures one of the 
approximations most commonly used in the electronic 
theories, the quasi-particle (QP) approximationpii, is se- 
riously questioned by the effect of the EP coupling. In- 
deed in most electronic systems characterized by a mod- 
erate internal correlation, the electrons are believed to 
occupy well defined energy levels characterized by a pre- 
cise energy, width and wave-function. The QP picture 
pictorially represents the effect of the correlation on these 
states as an electron-hole (in the case of electron-electron 
coupling) or electron-phonon (in the case of the EP cou- 
pling) pairs cloud which renormalizes the energy and the 
width of the electronic level, also reducing its effective 
electronic charge. The breakdown of the QP picture 
caused by the EP coupling has been already predicted 
in the case of superconductors by Scalapino et alJ^ and 



in complex metallic surfaces by Eiguren et alJ^. More 
recently we have sliowi*i^ a strong renormalization of the 
electronic properties of diamond and irans-polyacetylene 
caused by the EP coupling in the zero temperature limit. 

In this paper wc will extend our previous worki^, 
by providing more methodological and technical details 
of the dynamical theory we have previously used. We 
will also apply the same method to another polymer, 
polyethylene, finding a severe breakdown of the QP pic- 
ture. The analysis of the polyethylene results will con- 
firm and strengthen the general conclusions that we drew 
regarding the enormous impact of the electron-phonon 
coupling in carbon based nano-structures. 

In sections [IT] and IIIII we will review the derivation of 
the fully frequency dependent self-energy by using the 
many-body perturbation theory. The HAC theory will 
be, then, found as a static and adiabatic limit of the 
dynamical theory. In section IIVI and |V] wc will discuss 
how the structures appearing in the spectral functions of 
irans-polyacetylene and polyethylene rule out the basic 
assumptions of the HAC approach imposing the use of a 
fully dynamical theory. In section lVll we will show how 
the problem can be mapped in the solution of a fictitious 
Hamiltonian that makes possible to define the polaronic 
states as complex electron-phonon packets. Finally, in 
the conclusions, wc will point out as these results rep- 
resent an important step forward in the simulation of 
nanostructures, with a wealth of possible implications in 
the development of more refined theories for the elec- 
tronic and atomic dynamics. 



II. A DYNAMICAL APPROACH TO THE 
ELECTRON-PHONON PROBLEM 

We start from the generic form of the total Hamilto- 
nian of the system that wc divide in electronic (Hei), 
atomic (Hat) and electron-atom part [Hei-at)'- 



H — Hel + Hat + Hpj-at- 



(1) 



The Hamiltonian H admits both electronic and vibra- 
tional states that are coupled by Hei-at- In this work 
Density Functional Theory (DFT)i^ is used to calculate 

the eigenstates of H, where we use the notation O to indi- 
cate a quantity or an operator that is evaluated with the 
atoms frozen in their equilibrium crystallographic posi- 
tions. Similarly the vibrational states of the Hamiltonian 
H are described, fully ab-initio, by using the well-known 
extension of DFT, the Density Functional Perturbation 
Theory (DFPT>i^>ii. In DFPT the electronic correla- 
tions are embodied in a self-consistent mean potential 
Vscf representing the total electronic potential which de- 
pends on the atomic positions R/s = R/ -|- Ts : 



Hel-at 



I 



crystal 



drpir)Vsef[{R}]ir). 



(2) 



In the definition of R/^, / and s label the lattice cell (at 
position R/) and the atoms in the cell (at position r^), 
respectively. In Eq.([2]) p is the electron density operator. 

The aspect we are interested in this paper is how to 
properly include the modifications of the electronic lev- 
els induced by the atomic vibrations. In particular, by 
assuming the harmonic approximation for the phonons, 
we will develop a dynamical theory of the electronic dy- 
namics. To this end we follow a purely diagrammatic 
approacbi^ to present a short but accurate review of the 
derivation of the Fani^ self-energy and of the much less 
known Debye- Waller (DW) corrections^. 

If we now consider a configuration of lattice displace- 
ments uis, H can be expressed as a Taylor expansion 
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Isa^J s' (3 
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UlsaUja'P, (3) 



where a and j3 are the Cartesian coordinates. 

The link with the perturbative expansion is readily 
done by transforming Eq. ([3|) from the space of the lat- 
tice displacements to the space of the canonical lattice 
vibrations by means of the identity^i; 



iiisc. = Y. C^NgMsUJ^xr'/^ Ca (qA|s) e:^<^'+^^ 

qA 
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DqA 



(4) 



where Nq is the number of cells (or, equivalently the num- 
ber of q-points) used in the simulation and Mg is the 
atomic mass of the s atom in the unit cell. ^^ (q-^k) is 

the phonon polarization vector and ^LqA ^^'^ ^q-** ^^'^ ^^^ 
bosonic creation and annihilation operators. 
By inserting Eq. ([4]) into Eq. ([3]) we get 
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(5) 



In Eq. ([S|) we have introduced the first-order (gJJ/„i(.) and 

the second-order (A^,^^ ) electron-phonon matrix el- 
ements which will be shortly defined. To this purpose 
we rewrite Vscf making explicit its dependence on the 
atomic positions: 



Vscf[{R}]ir) = J2^scfii 



Ris). 



(6) 



From Eq. ^ it follows that the second order 
derivatives in the atomic positions are diagonal, 

^^-f^Vscf [{R}] (r) ex SijSss'^Vscf [{R}] (r). 



By using Eq. ([6]) the summation on R/ appearing in 
Eq. ([3]) leads to the momentum conservation both in the 
first and second order terms. At the first order this leads 
to the definition of the electron-phonon matrix elements 
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ynn'k 



^^^ -^(2Af,c.qA)-^/'e*'i- 



Ks) 



dV^ ' (r) 
X (nk| ;^^ \n'k - q)^^ (qA|.) , (7) 

We have also used the short form Rsa = Risa\j=o- A 
similar derivation can be followed to derive the 2"^^ order 
term 




(a) 2"'' order S^"" 




(b) 2"^* order E^ 



A qA.q', 



lEE 



C(qAk)C/3(q'A» 



X (nk| 



c'i^.^ai?. 
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'k-q-q'). (8) 



This second-order term is, in general, neglected as it is 
assumed to be small compared to the first-order term. 
Although this is correct at the level of the Hamiltonian 
it is not true anymore even at the lowest order of pertur- 
bation theory. 

Indeed the different terms in the Taylor expansion of 
the Hamiltonian defined in Eq. (O induce a wealth of di- 
agrams of increasing complexity and order. If we restrict 
to the lowest non vanishing order we have two diagrams: 
the Fan^ and the DW. These are presented by diagrams 
(a) and (6) in Fig. ([T]). In the same Figure two fourth 
order (in the displacements) diagrams, (c) and (d), are 
also showed. They are of the same order and, as the 
Fan and DW diagrams, they result from the perturba- 
tive treatment of the first order and second order terms 
in Eq. ^. 

The actual calculation of the Fan and DW diagrams is 
straightforward. The Fan's diagram is similar to the one 
generated by the electronic correlation in the so-called 
GW approximation!^, where the screened electronic in- 
teraction is replaced by a phonon propagator of wave vec- 
tor q and branch Ai^. Applying the finite temperature 
diagrammatic rules it is possible to define the Fan self- 
energy operator S^^" (*'^ii^); recovering the expression 
originally evaluated by Fan^^: 



X^I< an 



(-^^r) = -^]^EEi5.'?«'.i'x 



-|-OC' 



qA n' 

-.(0) /„■, , ^ ^(0) 



E Kxi'^^)G^ni-^(^u,-iu,), (9) 



j = -oo 



where /3 = -^ [k is the Boltzmann constant) and T 
is the temperature of the phonon bath. By using the 
standard definitions of the electronic and the phononic 



Green's functions: G 
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(where /i is the chemical potential), 
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order E 




(d) 4"^ order E 



FIG. 1. The self-energy diagrams corresponding to the first 

and second order terms in the Taylor expansion oi H — H (see 
Eq. l[3]l) treated at different orders of the perturbativc expansion. 
For example the well-known Fan self-energy is formally obtained 
as a 2"'^ order expansion of the first order term in Eq. ^. However 
the second term of Eq. l[3} , treated at first order gives the 2"'' order 
j^DW ^ that is of the same order of the Fan term and, consequently, 
cannot be neglected. The diagram (c) is obtained as a 4*'' order 
expansion of the first order in Eq. {Sj while (d) comes from the 
second term of Eq. |(3} treated at the 2*'' order. 



( 



{iojj — iOqx) — {i^j + ^qA) ) , and summing over the 
Matsubara frequencies, we get the final expression for the 
Fan self-energy 



i],^r(«^,T) = E 



i^nn'k 



n'qA ^ 



AfqA (T) + 1 - fn'k-d ^ 

iuj - e„'k-q - WqA - iO+ 

, A^qA (T) + /„.k_q 

iu - £„'k-q + UJqX - «0^ 



(10) 



where N^x {T) is the Bose function distribution of the 
phonon mode (q. A) at temperature T. 

A similar expression can be derived for the fre- 
quency independent DW self-energy S^j^ . This 
term comes from the equal time contractions of the 

^qAJ f^-q'A' +K'X') operators: 



iU 



{{bl 



qA 



'qAj [b-q'X'+K'y)) = 

(5_q,q-(5A,A' [Aq^A (T) + AqA (T) + 1] . (11) 



The corresponding diagram (6) in Fig.[l] can be easily 
found to be 

^n^ iT)-^T. ^™ir''' (2iVqA (T) + 1) . (12) 



qA 



Both the Fan and DW self-energy have been already de- 
rived previously in the framework of the Heine-Allen- 
Cardona (HAC) theorjii^i^. The HAC approach is based 
on the static Rayleigh-Schrodinger perturbation theory. 
More precisely the uisa are used as scalar variables on 
which a static perturbation theory is applied. As we will 
mention in Sec lIIIi the second-order derivatives appear- 
ing in the definition of the DW term can be rewritten 
in terms of the one-order derivatives by imposing the 
translational invariance of the correction to the electronic 
levels. 

On the other hand the Frohlich and the Holstein 
Hamiltonians usually neglect the DW term (diagram (6) 
in Fig.[T]), even if it is of the same order of the Fan term. 

The many body formulation represents the dynamical 
extension of the HAC approach, that is recovered from 
Eq. (|10p by using uj w £„k (the on-the-mass-shcll(OMS) 
limit) and |£„k — £n'k-q\ ^ ^qX (the adiabatic limit) and 
by considering only the real part of the self-energy. It 
turns out, therefore, that in the HAC approach the tem- 
perature dependent change in the single-particle energies 
is given by 



Ae^k^^(r) 



^DW 



(T) 



^r-kl 2iVq,(T)-l-l 



n'qA 



Na 



£^nk S^n'k- 



(13) 



The soundness of the HAC approach is then, from a 
MBPT perspective, connected to the validity of the on- 
the~mass-shell and of the adiabatic approximations. We 
will prove in the section lTVl that these approximations are 
not always well motivated. Dynamical and non-adiabatic 
corrections can be huge, such to invalidate the applica- 
bility of the HAC approach. 



III. SECOND ORDER DERIVATIVES OF Vscf 

AND THE ACTUAL CALCULATION OF THE 

DEBYE WALLER SELF-ENERGY 

The general expression, Eq. ([5]) and Eq. ([8|) for the 
perturbed Hamiltonian requires the knowledge of the 
second-order gradients of the self-consistent potential 



n'p,nk 



(npl 



d'v^:} 



(r) 



dRs^dR 



sp 



rik). 



(14) 



where p here replaces the k -I- q -I- q' vector appearing in 
Eq. ([5|) and Eq. (j5]) . These terms are extremely cumber- 
some to calculate and the task becomes easily prohibitive 
when higher orders are included. Their evaluation is. 



however, crucial because the k = p case is needed to cal- 
culate the lowest-order DW self-energy, while the finite 
momenta matrix element defines higher order diagrams 
(like diagram (d) in Fig. ((!])). 

In the case of the simpler S^i^ (diagram [h) in Fig. ([1])) 
we know, from Eq. ([8]) and Eq. (fT4| that 



A:!::r' = jEE(2^^^«^'iA)"' 



a, [3 



C(qAk)e/^(-qA|g) 

2A/sWqA 



A 



tk.nk' 



(15) 



In order to evaluate the A!,?f „. factors the HAC the- 
ory uses the fact that if all atoms were shifted by the 
same amount all physical quantities should not change. 
In other terms, being the Ae^j^*^ (T) an explicit func- 
tional of the atomic positions (Eq. (JH)) that are treated 
classically, it is possible to impose the following transla- 
tional invariance condition 



Ae^k^^ 



{uisa}] (T) 



Aefk^^ 



[{uisa+da,}]{T). (16) 



From this condition it follows thatii^i^ in order to cal- 
culate A^^j^'* that defines the DW self-energy (see 
Eq. p^ *) only the matrix clement A^'^'^j^. is needed. This 
can be rewritten as 



AT,. 

nk.nk 



E 



^nk ^?i'k 



Z^("k|^j^^ |7ik) I {nk\ ^^ Ink) 



dRs'a 



OR 



s/3 



-(nk| — ^ ink) > (nk| — ^ \nK) 



dRsa 



dRs't 



(17) 



The condition given by Eq. (|16|) is, however, intrinsically 
ill-defined in the diagrammatic approach: the correc- 
tion to the energy levels is a quantity obtained in fact, 
from the self-energy operator that, in turns, can be de- 
fined only when the displacement operators are quantized 
and a second quantized form of the Hamiltonian change 
(Eq. (O) is introduced. 

This inconsistency can be, indeed, cured in a fully 
MBPT framework^S but it requires to introduce the con- 
stant displacement vector d^ as an operator. This leads 
to the definition of new kinds of diagrams that will de- 
pend on powers of da ■ By imposing that diagrams of the 
same order cancel each other it is possible to obtain a 
general expression for the matrix element A^?^ ^^j^. 

As discussed by X. Gonze^^, the local dependence on 
the atomic positions in Eq. ^ assumes that the elec- 
tronic screening of the ionic potential, that defines Vscf, 
depends only smoothly on R/s. By taking fully into 



account this intrinsic dependence on the atomic posi- 
tions a correction to the Debye- Waller term, named non- 
diagonal Debye-Waller correction, can be defined. This 
correction has been reported to be important for isolated 
molecules and atoms^. Its effect in solids and, more gen- 
erally, in extended systems is expected to be weakened 
by the efficient screening properties. 



constant, if the bare energy e„k is far from a pole of self- 
energy, then S^^" i^'"^) ^^^ be Taylor expanded, up to 
the first order, around Snk'- 

En,. (T) = e„k + S,T (T) + S^r (£nk, T) + 



dKr 



{^,T) 



duj 



(i^nk (T) - e„k) . (21) 



IV. DYNAMICAL SELF-ENERGY EFFECTS 

BEYOND THE QUASI PARTICLE 

APPROXIMATION 

In section|lT]we showed that the HAC theory represents 
the static and adiabatic limit of the dynamical electron- 
phonon self-energy. The more suitable are the conditions 
of validity of the on-the-mass-shell and of the adiabatic 
approximations, the sounder is the applicability of the 
HAC approach from a MBPT perspective. We want to 
prove in this section that these approximations are not 
always well motivated and dynamical and non-adiabatic 
corrections to the HAC approach cannot be neglected a 
priori. 

The HAC approach grounds on the concept of a well 
defined QP state: the charge carriers are assumed to be 
concentrated on electronic levels, being characterized by 
a well defined energy and wave-function. The QP con- 
cept can be firmly introduced in a many-body Green's 
function theory— where the definition embodies, at the 
same time, its limitations, as it will be clear in the fol- 
lowing. 

The fully interacting propagator can be written, for 
real energies w in terms of the self-energy, (Eq. ©) as 



G„i.{u),T) 



1 



f nk 



^Fan 
^rik 



i^,T) 



^nk 



(T) 



(18) 



The rotation from the imaginary to the real axis has been 
easily performed by a Wick rotatioi*2i as the energy de- 
pendence of the Fan self-energy is explicit. The single 
particle excitations are then the complex poles of Eq. p^ 



-z5[E^r(i?nk(T),r)]=0. (19) 



As it is clear from Eq. p9|) a genuine QP state should 
have a zero line-width, that is a zero imaginary part of 
the self-energy. In practice this is never completely true. 
Nevertheless, when the frequency dependence of the self- 
energy is smooth Eq. (|18p can be rewritten by using two 
simple and intuitive approximations: the OMS and the 
QP approximation. In the specific case of a constant and 
real self-energy (as in the HAC case) one can introduce 
the OMS approximation where the solution of Eq. P^ is 
given by 



Enu (T) = £„k + S,T (T) 



Y^_^ an 
^nk 



£7ik,T) . 



(20) 



We notice, from Eqs. dH]) and ^, that Enw{T) is 
complex. Even in the case where the self-energy is not 



Eq. (pij) corresponds to the QP approximation. The bare 
energy is then renormalized because of the virtual scat- 
terings which are described by the real part of the self- 
energy. This renormalization is easily described by the 
solution of Eq. (PT|) : 

En]i (T) = £„k + 

+ Z„k (T) [E^k"" (^»k, T) + J:^^ (T)] , (22) 



with Z„k (T) 



asS"('^,r) 



From Eq. 



-1 

the renor- 

UJ—Sn 

it is evident that 



malization factor. From Eq. (|22] 
Enk is complex and its imaginary part r„k (T) ~ 
'^[EnkiT)]^ the QP line- width, is proportional to 
the S[Z„k(T)E^r(£nk,T)]. A smaU r„k (T) indi- 
cates a stable QP, that slowly decays because of the 
real scatterings with the other particles and with the 
phonon modes. By assuming the QP approximation 
to be valid Eq. (fTS)) can be re-written as G„k {<^, T) = 

ZnkiT){LU-Enk{T)y\ 

Angle Resolved Photoemission Spectroscopy (ARPES) 
provides a definitive tool to verify if the QP approxima- 
tion is accurate. If it existed a true QP should appear as a 
peak in the photoemission spectra with a Lorentzian line- 
shape. Indeed one finds that the spectral function (SF) 
Anw{oJ,T) = TT-i I 3[G„k(w,T)] I is given, in the QP 
approximation and in the simple case of a purely real 
^nk, by 



4f (--7^) 



Znk (T) |r„k (T) I 



n {oj-^[Enu{T)]f + rl^{T) 



(23) 



In this case the QP energy and width give the peak po- 
sition and the spectral peak width. It is worth noticing 
that a complex value of ZnW (T) would cause the spectral 
function to have an asymmetric lineshape. 

The SF gives a physical interpretation and a clear val- 
idation of the QP approximation. Indeed A^^' (w,T) is 
a probability function to find an electron in the state nk. 
with energy w and the total electronic charge associated 
to the QP state is Zni.{T), that corresponds to the in- 
tegral of ^„k i^^T). The renormalization factor repre- 
sents, therefore, the QP charge. When Znw {T) = 1 and 
3 [Snk(enk,T)] — > the SF reduces to a delta function, 
the SF of a particle with energy K [Enk (T)]. 

It is clear that a direct comparison of ^Jfk ('^i ^) with 
the true SF corresponding to a given self-energy or with 
the ARPES lineshape provides the ultimate validation of 
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FIG. 2. Trans— polyacetylene. Spectral function (upper frame) and 
self— energy (lower frame) corresponding to the state | n = 1 , k = 
r) . In the self— energy frame both the real (solid line) and imaginary 
(dashed line) parts are showed. The three arrows represent the 

bare electronic energy e^k and the two solutions (^„k'^nk) °^ 
Eq. I|19| l. The thin solid straight line represents instead the function 
t^ ~ ^nk ~ ^nk^' ^^ *^^ imaginary part of the self— energy shows a 
clear, intense and wide peak at around — 16.5eV (i.e. very close to 
Ejik) the real-part is dominated by a rapid oscillation that cannot 
be captured at all by the linearization of the energy dependence 
and causes the appearance of two solutions of Eq. (|19ll . 



the QP picture. A paradigmatic example that well ex- 
plains the basic mechanism for the breakdown of the QP 
picture is given in Fig. [2] where the zero temperature SF 
for the I n = l,k = F) state of trans-polyacetylcnc is 
showed in the upper frame. It is clear that the SF of this 
state is far from being well represented by a Lorentzian 
lineshape. Indeed it is evident the appearance of two 
peaks at energies E^^^^ and E\2^. These two peaks arc 
the signature of a breakdown of the QP picture, because 
the existence condition of only one pole collecting most 
of the weight is not satisfied. We will come back on the 
physical interpretation of such structures in the next sec- 
tion. Nevertheless the origin of these two peaks is evident 
if we analyze the energy dependence of the imaginary and 
real parts of the corresponding self-energy, showed in the 
lower frame of the same figure. In this specific case the 
solution of Eq. (fTO)) admits two roots as a consequence 
of the rapid oscillation of 3? [E^^" (w,T = 0)] around 
— 16.5eV (the bare electronic energy of this state). This 
oscillation is, in turn, induced by an intense peak ap- 



pearing in the 3 






^°" (w, T = 0)] . We deduce that in 
this case a naive application of the QP approximation in 
form of a linearization of S^^" (w, T) (Eq. (PT|) ) may pro- 
duce a non physical energy dependence of the self-energy. 
As a consequence this leads to meaningless (negative or 
enormously large) values of .^„k- 

In contrast to the case of purely electronic self- 



energies the QP approximation is known to lead to a 
too rough description of the electron-phonon spectral 
function for low-energy electrons in the homogeneous 
electron gas (jellium). As discussed by Engelsberg and 
Schrieffer— this failure, although not as dramatic as the 
one found in the present case, is linked to the mixing of 
electronic and phononic excitations. More precisely the 
authors identify three kind of excitations that appear as 
poles of the Green's function when the energy of the elec- 
tronic level is increased well above the Debye energy and 
the strength of the electron-phonon coupling is also in- 
creased. One is a purely QP state where the electron is 
dressed by a phonon cloud. The others lie in the contin- 
uum of electron-phonon pairs composed by clothed elec- 
trons and clothed phonons being excited, with a constant 
momentum sum. 



V. BREAKDOWN OF THE QUASI PARTICLE 

APPROXIMATION: THE CASE OF 

TRANS-POLYACETYLENE AND 

POLYETHYLENE 

In the previous section we showed that the structures 
appearing in the SF of frans-polyacetylene rule out any 
description of the coupled electron-phonon system in 
terms of QPs. More importantly the rich structure of 
peaks appearing in the SF described in Fig. [2] is not a for- 
tuitous case. It is actually a general trend both in trans- 
polyacetylene and in polyethylene. Indeed, in Figs. [3] and 
21 the bare electronic band structure and the correspond- 
ing SFs for a fixed k vector are shown in the upper frame. 
The position of the k vector in the Brillouin Zone (BZ) 
is represented by an horizontal line in the lower frame 
where the the valence bands of the two polymers are also 
reported. 

In the upper frame of Figs.[3]and|3]there is also a sketch 
of the atomic structure of the polymers that helps to un- 
derstand the main key differences among them. Both 
systems are linear polymers. Trans-polyacetylene is a 
conjugated polymer where each carbon atom forms four 
nearest-neighbour bonds. Three of the four carbon va- 
lence electrons are in sp^ hybridized orbitals and two of 
the cr-type bonds connect neighbour carbons along the 
one-dimensional backbone, while the third forms a bond 
with the hydrogen side group. Polyethylene , instead, 
is a tr-bonded, non-conjugated polymer. The atoms in 
the unit cell do not lie on a single plane like in trans- 
poly acetylene. The C atoms are sp"^ hybridized and, as 
in the polyethylene each C atom has four bonds. 

From the upper panels of Figs.|3]and|4]it is evident that 
the EP interaction dramatically affects the spectral func- 
tions. The most striking aspect is that the SFs exhibit a 
multiplicity of structures. Although in some cases a sin- 
gle and strong peak can be observed the general trend is a 
very complex ensemble of peaks that makes impossible to 
apply the QP approximation. We will give a more formal 
and mathematical description of the internal structure of 
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FIG. 3. (color on line). The SFs of the last four occupied states 
in irans-polyacctylene (upper frame) are shown together with the 
DFT-LDA bands (lower frame). The horizontal line in the bands 
frame represents the position of the k-point and the energy range 
along which the SFs are displayed. The vertical arrows in the 
upper frame represents the energy position of the unperturbed DFT 
levels. Since these states correspond to in plane orbitals they are 
strongly affected by the in plane atomic vibrations. The result is 
that the bare electronic levels are split in several polaronic states. 



these peaks in the next section. Here we would Uke to 
underline some of their general aspects. 

A remarkable aspect is that some SFs are so largely 
structured that they span a large energy range, even 3 eV 
(see the 6*'* band of polyethylene, Fig.|4]). If they span 
a so large energy range, the SFs end up with overlap- 
ping each other in some cases (like the 4*'' and 5*'' band 
SFs in trans-polyacetylene) . The crucial and straight- 
forward consequence is that it turns out difficult to as- 
sociate a single and well defined energy to the electron 
and, more importantly, different bands will energetically 
merge pointing to a non trivial mixing of the electronic 
states. 

When a SF covers a large energy range, one peak may 
be distant in energy from the others more than the Debye 
energy (0.4 eV in irans-polyacetylene and polyethylene). 
For this reason each peak can not be simply interpreted 
as a main QP peak plus a phonon replica. This point will 
be further discussed in Appendix A by using a two band 
model. Nevertheless it is reasonable to speculate that the 
formation of more than one peak suggests to reformulate 
the problem in a different framework where bare elec- 
trons are mixed with phonons, and not simply screened 
by phonons. Each peak appearing in the SF would be 
then identified by a mixed electron-phonon states. Since 
the many body framework is not suitable to add infor- 
mation about the composition of the new mixed states, 
we will reformulate the problem in the next section by 
mapping the problem into and Hamiltonian representa- 
tion. 



FIG. 4. (color on line). Like in Fig.ljSJI in the polyethylene case. 







FIG. 5. (color on line). Two dimensional plot of the spectral func- 
tions AZ„i^{uj) for polyethylene in the last four occupied bands re- 
gion. The DFT and the polaronic bands are opportunely labeled. 
In general the electronic levels acquire a large energy indetermina- 
tion if compared to the DFT bands represented by solid black lines. 
The EP interaction moves up of about 300 meV the last two occu- 
pied bands leading to an increase of the band width. Moreover the 
6"" band near the X point shows a large energy indetermination 
that makes it almost disappearing. 



A global view of the effect of the ZPM on the electronic 
structure of polyethylene is given in Fig. [S] in the energy 
range of the last 4 occupied valence bands. The DFT 
electronic bands are drawn as a reference of the electronic 
band structure before switching on the EP interaction. 

By defining AZ„k (w) = Ank (w) Aw the probabil- 
ity to find an electron | nk) in the small energy range 
Aoj = 50 meV, we made a bidimensional representation 



of the probability amplitude. This is showed in Fig. [5] 
by using a colored scale that goes from white (the less 
intense peak), to black (the most intense one). As a con- 
sequence this picture gathers all the information about 
the energy range covered by the SFs and the intensity 
of all peaks. In particular we observe that the 6*'' band 
of polyethylene moves up close to F-point and then the 
electron completely disappears. 

The resulting zero point renormalization of the gap at 
F point of polyethylene is 280meV, larger than the trans- 
polyacetylene casei^. Such a difference is ascribed to the 
peculiar shape of the irans-polyacetylene orbital at the 
X point whose tt character corresponds to states per- 
pendicular to the polymer axis. Thus they feel less the 
effect of the in-plane vibrations. On the other hand for 
polyethylene at F, the electrons arc localized along the 
C — C bond, where the zero point motion effect of the 
electronic gap is sizable. For what concerns the deeper 
states far from the gap, the effect of the electron-phonon 
coupling is equally strong. In fact as they are in plane 
orbitals they are directly affected by in plane atomic vi- 
brations. 

We also observe that each band has a different energy 
width which evolves in different manners moving from 
T io X. An increasing of the bandwidth is normally 
associated to a consequent increase of the delocalization 
of the orbitals. This fact can be used link the effect of the 
EP coupling to a increased electronic mobility mediated 
by the polaronic states. 

In the next section we will go beyond this picture by 
introducing a general framework to link the poles of the 
electron-phonon Green's functions to coupled packets of 
electron-phonon pairs. 



VI. INTERNAL STRUCTURE OF THE 

POLARONIC STATES VIA AN HAMILTONIAN 

REPRESENTATION 

In order to gain more insight into the complex struc- 
tures that appear in the SFs of both irans-polyacetylene 
and polyethylene we propose, in this section, a mapping 
of the Many-Body problem into an equivalent Hamilto- 
nian representation. In this representation the poles of 
the SF will appear as eigenvalues of a fictitious el-ph 
Hamiltonian. When this is solved in a specific restricted 
sub-space of the entire Fock space, it will reproduce the 
same SF obtained by solving the Dyson equation within 
the Fan and Debye-Waller approximations for the self- 
energy. 

In order to show this we start by rewriting the SF by 
using the well-known Lehmann representationi^: 



A„k (c., T) = 5^ |(vI/o I 4^ I /k (T)) \^5 {u - Ej^ (T)) , 

(24) 



/k 



and real poles of the GF. As our initial Hamiltonian, 
Eq. (dl , is composed of electrons and phonons the states 
Ik (T)) live an extended Fock space composed of elec- 
trons and phonons. 

In the QP approximation the distribution of peaks ap- 
pearing in Eq. (j24p is approximated with a Lorentzian 
distribution centered at the QP energy with a width 
equal to the QP line-width. Thus, Eq. (pi)) already 
underlines that the origin of the multiple poles in the 
SFs shown in Fig.|3] and Fig.|3] is connected to the exis- 
tence of more than one intense state | /k) belonging to 
the same state | nk). 

Now we assume that Eq. (IMl) remains valid also when 

the exact self-energy is approximated by the Fan and 

Debye-Waller terms. And, in order to link the states 

/k (T)) to an Hamiltonian problem we start rewriting 

Eq. ([1]) in second quantization: 



H — Hel + Hpii + Hel- 



el—ph: 



(25) 



where Hei is the electronic Hamiltonian, Hph is the in- 
dependent phonons Hamiltonian and H^i-ph is the EP 
interaction Hamiltonian. The last three terms, written 
in the second quantization, read 



Hei = ^ e«kc|,kCnk, (26) 

nk 

Hph^Y.^cx{h{j^c^x + \), m) 

\- E ffn„'kc|;kC„'k-q(&^_q,+6qA). m) 



H el-ph 



N, 



n.n .k 
q,A 



where | Jk (T)) arc the true cigcnstates of the system 
with energy Eji^ (T) which, in turns, represent the true 



In Eq. ([26|) £„k is a single particle energy that we will 
shortly define, c^j^ is the creation and c„'k-q is the anni- 
hilation electronic operators, b^ and bq\ are the creation 
and annihilation operators for phonons with energy ujq\ 
and wave vector q. 9^'^/]^ are the EP coupling matrix 
elements (see Eq. (O). 

We want now to use Eqs. (|^5|) to calculate the states 
|/k(T)). To this end we note, from the (a) frame of 
FiglU that the Fan self-energy makes an initial state 
|nk) to scatter with a phonon state |qA), with population 
Nq\ ± 1, in a final state |n'k — q). Only one phonon is 
exchanged and in the self-energy loop the intermediate 
states are the composite pairs | n'k — q)® | Nqx^^) with 
energy e„'k-q±WqA. This energy, indeed, appears in the 
denominator of Eq. ([TU| . 

Physically this means that, if we introduce the gen- 
eral state product of an electronic and phononic part 

n'k — q)(8) | Nq\ ± 1), the intermediate states of the 
self-energy arc all possible combinations with different q 
and A. It follows that we can guess, at a given tempera- 



ture, 



\IkiT))=Y,Ai^iT)\nk) + 

n 

+ E ^"'k-q {T) I n'k - q)® I iVqA (T) ± 1). (27) 

n'qA 

The coefficients A^ and B^^ can be found by diago- 
nalizing Eq. (|25|) in the space of electron-phonon states 
spanned by the definition given in Eq. (P7)) . More pre- 
cisely, in order to expand the matrix form of Eq. (|25p we 
notice that, due to Eq. ((27| the basis set will be composed 
of the following elements 

Ink) liVqA), |nk-q) | A^qA ± 1). (28) 

At zero temperature the basis set is reduced to 

|nk)|Op^), |nk-q)|lqA), (29) 

and it reflects the fact that at T = QK there are no 
phonons in the ground state. In this restricted basis the 
Hamiltonian reads 



H = 



£^nk 



^nn' 



rjep 
nn'qA 



^-"rm'qA/ 



Snk-qSnn' + <^ci\Sqq'S\X' 



(30) 



The fact that Hei-ph is Hermitian makes also H Hcr- 
mitian. The equivalence of Eq. (P^ with the spectral 
function calculated in Sec lIVI is obtained by imposing 
£nk = Erik + ^Sc^' This equivalence is proved analyt- 
ically, in the zero temperature limit, in Appendix 1X1 for 
a simple two-levels model. 

Prom Eq. (j30|) we also note that the number of states 
I /k(T)) is equal to the dimension of the matrix and 
it is obtained by multiplying the number of electronic 
bands times the number of q vectors times the number 
of phononic branches, A. As a consequence the number of 
states I /k (T)) is larger than that of Hei- This confirms 
the fact that, in Eq. (|24p . the states |/k (T)) form a con- 
tinuum that, in the QP approximation, dresses the bare 
electronic state. This dressing is, de-facto, represented 
by a cloud of mixed clectron-phonon states surrounding 
the QP energy with a Lorcntzian distribution. 

The Hamiltonian H (Eq. ((30l) l is diagonalized for 
frans-polyacetylene including 30 electronic bands, 10 q- 
vectors and 12 phonon branches. Since the Hamiltonian 
is written on a complete basis, the coefficients A^j, and 

i/A 
n'k— q 



Bi}, satisfy the following condition 



E 



At 



E 

n'qA 



^n'k-q 



1, 



(31) 
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FIG. 6. Trans-polyacetylenc. The SF of the state | n = 4, k = 
0.2(— 2^, 0, 0)) is decomposed in polaronic states, each labeled by 
I /k). Several structures appear thus ruling out the QP approxi- 
mation. 
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FIG. 7. Trans— polyacetylene. The projection of the polaronic 
states over the corresponding pure electronic state, Z^^j^ is shown. 
The dashed line represents Z^^ of a pure electron state as a refer- 
ence value. 



which ensures that the spectral function A„k(a;,r) is 
correctly normalized to 1 when integrated over all fre- 
quencies. Once the eigenstates | /k) and the eigenvalues 
i?/k are known, the SF can be calculated according to 
Eq. ((M)) and all the peaks appearing in the SPs of state 
nk), are unambiguously labeled with a particular | /k) 
state, having | nk) as the pure electronic component. Let 
us consider the | n = 4,k = 0.2( — ,0,0)) state of trans- 
poly acetylene as an example. In Fig.|6]it is shown the cor- 
responding zero temperature SF. The poles and the cor- 
responding residuals are indicated in Fig. [HI by bars with 
different heights. The residuals are given by | A^^ p, 
that is the probability to find the polaronic state in the 
pure electronic | nk) state. This reminds the physical 
meaning of the Znk factors, and we use this similarity to 



define Z.^^ =| A^^ |^. Nevertheless from Eq. ^7^ it is ev- 
ident that the smaller the Zj^^^ is, the less the polaronic 
state can be assimilated to an electron. It means that 
the mixed EP contribution in Eq. (P7)) indeed weights the 



Ai 



can even be as 



most. In the case of Fig. |6] the 
small as 0.2. 

These small values of Z^j^ represent a general trend. In 
Fig.[7]Z^jj. is plotted as a function of the polaronic eigen- 
values. It can be noted that only few polaronic states 
have Z,{j^ ~ 1. Most of all are below 0.5, instead. It 
means that the mixed EP part of the eigenstate, shown 
in Eq. ()27p , plays a dominant role. 

A small value of Z^j, points to a non trivial physical 
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Trans-polyacetylene 


polyethylene 




C 
a.u. 


H 

a.u. 


C 
a.u. 


H 

a.u. 


X 


0.18 


0.55 


0.1 


0.32 


y 


0.13 


0.36 


0.07 


0.21 


z 


0.11 


0.56 


0.07 


0.34 



TABLE I. Atomic amplitudes obtained by evaluating the matrix 
elements of operator \/u|jj.^. Because of its smaller mass the hy- 
drogen quantum atomic size is three times larger than the carbon 
atom one. Nevertheless the constraint imposed by the different 
geometries makes the deviation of the two polymers appreciably 
different. 



property of the polaronic states. When Z^^^ ^> it is 
evident that only the mixed terms in the sum, where the 
electrons and the phonons appear together, are non zero. 
This means that the electrons cannot move in the sys- 
tem alone, even if dressed by an electron-phonon cloud, 
but need to build up true bound clcctron-phonon states. 
This is a clear fingerprint of the breakdown of the QP 
approximation. 



The generic definition of the polaronic state, Eq. ([27| . 
allows to calculate the mean value of any observable that 
lives in the mixed electron-phonon space. For example 
we can evaluate the matrix elements of the atomic inde- 
termination operator as follows 



"a/ks = (-^k I U^, I /k) = 
1 



E 

qA 



2NqUj^xMs 



Ea (qA/s) el (qA/s) x 



EM"kP+3El^"'k-q 



(32) 



By using these m^/^s ^^ "-^^^ associate an average quan- 
tum size to the atoms. The values for the C and H 
species, calculated by Eq. ([32]) . are shown in Tab. HI These 
values point to the fact that the atoms acquire an indcter- 
mination larger along the polymer axis. Since H is lighter 
than C the atomic quantum size is larger. The different 
constraint created by the geometries is the cause of the 



different ^/w-^ 



:7k; 



values between the two polymers. 



The values of the atomic indctermination suggest that 
electrons and phonons exert a cooperative effect on each 
other. The charge density spreads all along the polymer, 
while the atoms squeeze along y and z directions, widen- 
ing along X. This cooperation can cause, for example, 
an enhancement of the mobility, opening therefore new 
perspectives for future investigations and applications of 
polymers. 



VII. CONCLUSIONS 

In this work wc have shown that the Hcine-Allen- 
Cardona approach suffers of some severe limitations when 
applied to predict the zero temperature energy correc- 
tion in low dimensional systems. We extensively de- 
scribe a fully dynamical extension of the Heine-Allen- 
Cardona approach to show that the zero point motion 
effect severely questions the reliability of the QP picture 
in irans-polyacetylene and polyethylene. 

The single particle spectral functions, indeed, exhibit 
multiple structures at T = /v . The formation of addi- 
tional structures caused by the strong electron-phonon 
interaction is interpreted in terms of composed electron- 
phonon states, what we call in this work polaronic states. 
These states are precisely defined by mapping the struc- 
tures of the Many-Body spectral functions into the solu- 
tion of an eigenvalue problem. 

Thanks to this important mapping the non perturba- 
tive nature of the polaronic states appears as a coherent 
superposition of clcctron-phonon pairs. And the coop- 
erative dynamics between electrons and atoms in these 
states rules out any description in terms of bare atoms 
and quasiparticlcs. 

The resulting coupled electronic and atomic dynam- 
ics pave the way for new investigations in polymers and 
more in general in low dimensional nanostructures. The 
cooperative dynamics of electrons and phonons in the po- 
laronic states can have potential physical implications, as 
for example, an enhancement of the electronic mobility. 

More generally the breakdown of the quasiparticle pic- 
ture imposes a critical analysis of the previous results 
obtained using purely electronic theories. 



Appendix A: A two— levels model to verify the 
Hamiltonian representation 



In Sec lVIl we have given physical arguments to sup- 
port the choice of a specific limited Fock space where the 
Hamiltonian problem is solved. These arguments were 
guided by the final goal of introducing an Hamiltonian 
representation that gives exactly the Green's functions 
corresponding to the Fan approximation for the self- 
energy. To better investigate and confirm this ansatz 
let us consider two levels of energies e^ = 0, -E coupled to 
a phonon of energy wq at T = A'. 

Eqs. (|25ll30p thus reduce to a simple expression for the 
Hamiltonian of this system 



H 



V e,c\c. 



- ^ob% + Y^ gc]c^{bl + bo). (Al) 



4=1,2 

i=l,2 



As discussed in ScclVIIwe consider for the finite basis the 
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following ansatz: 



|l)|0p^),|2)|0p,), 
|1)| V),|2)|lp,), 



(A2) 
(A3) 



where by | i) we mean an electron in the level i*''. The 
dimension of the Hamiltonian matrix is then given by 
multiplying 2 bands x 1 q point x 1 phonon branch, re- 
sulting in a 4 X 4 matrix 



U 



fo g \ 

E g 

g ujo 

\g E + Looj 



(A4) 



which can be diagonalized in two blocks, obtaining the 
following four energy levels 



Ei = 



Eo 



E?, 



E + uja + 


^{E + 


^0? 


+ V 




2 






E + uja- 


^J{E + 


^^Y 


+v 




2 






-E + WO + 


^{E- 


^0? 


+v 




2 







(A5) 



(A6) 



(A7) 



E, . ^ + -o-V(^--o)^+_V. (A8) 



The corresponding four eigenvectors are 



/4) = T^ 



-/V4 \ Wq — -B; 



,1,0 , 



(A9) 
(AlO) 
(All) 
(A12) 



where Ni are the normalization factors, with A^i = N2 
and Nz = Ni. 

These are the needed ingredients to calculate the 
Green's functions as matrix element of the resolvent 

Qi{uj) = {vac I ci — c+ I vac), (A13) 

(jj — H 

where | vac) is the vacuum of phonons and electrons. Ex- 
panding in eigenstates of the system, Eq. (JAISP becomes 

Gii^) = "^{vo-c I c, _ I Ij){Ij I c+ I vac) 



4 



oj-H 



51 1 ^"""^ I ^' I ^j) I^TTT' 



j=i 



oj — Ei 



(AM) 



In order to show that the SFs calculated from Eq. (|A14[) 
are equivalent to the ones obtained in the MB approach, 
the Green's function for the f** state is evaluated from 
Eq. (|A14p as follows 



Qiiuj) 



1 



El — E2 



'E, 



El 



OJ — El u) — E2 



(A15) 



On the other hand the Fan approximation for the self- 
energy, Eq. pH)) . in this test case reduces to 



Si^ an f \ A 

1 K^)^g 



iv(c^o) + 1-/2 



E- 



ujQ — zO+ 

iV(wo) + /2 



E + uja- tO+ 



(A16) 



At zero temperature the Bose occupation factors vanish. 
The Fermi occupation factor /2 is zero because the level 
is empty and Eq. (|A16p becomes 



GiH 



^ - g 



uj — E—ujQ 



iO+ 



(A17) 



The poles of Eq. ()A17j) are uj — Ei and uj = E2, defined 
by Eqs. (|A5HA6|) . The residues evaluated at each pole 
are shown as follows 



Lj = El Resi = 
UJ = E2 Res2 = 



E-y 



Ei- 

Ei 



E2 



El — E2 



(A18) 
(A19) 



The final expression for Qi in the many body approach 
is then 



Giioj) 



1 



E\ — Eo 



E, 



El 



oj-Ei- iO^ 



UJ-E2- iO-i 



that is equivalent to Eq. (|A15[) . 
From Eqs lA5f[A6l we notice that 



Ei^E2^ yJ{E + wo)2 + 4,92 



(A20) 



(A21) 



which is larger than ojq and it is also larger then the same 
energy difference when g — >■ 0. This clearly means that 
both electrons take part in the formation of the polaronic 
state thanks to the additional energy provided by the EP 
coupling. This also implies that each additional structure 
cannot be interpreted in energetic terms as simply an 
electron "plus" one phonon. 



Appendix B: Calculation Details 

The phonon modes and the electron-phonon matrix el- 
ements were calculated using a uniform grid of 10 x 1 x 1 
k~points. We used a plane-waves basis and norm con- 
serving pseudo-potentials— for the carbon and hydro- 
gen atoms. The exchange correlation potential has been 
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treated within the local density approximation. For the 
ground-state calculations we used the PWSCF code—. 
The Fan self-energy and the Debye-Waller contribution 
are calculated using a random grid of transferred mo- 
menta, using the yainbo code—. 

The numerical evaluation of Eq. ([T0| is a formidable 
task. Indeed the use of a fine sampling of the BZ is 
prohibitive. The reason is that such large grids of trans- 
ferred momenta are inevitably connected with the use of 
equally large grids of k-points. An alternative solution, 
that we used in the present calculations, is to fix a cer- 
tain k-points grid and to perform the integration of the 
BZ by using a random grid of points to perform the q 
summation in Eq. (jlOp . 

Moreover, in order to speed-up the convergence with 
the number of random points and to take in account 
the divergence at q — > of the |ff^'„kP matrix elements 
we divide the BZ in small spherical regions Rq centered 
around each q point. Eq. (|10[) can be then rewritten as 



Q n'\ V^^Q 



qA 

ynn'k 



Nqx iT) + l- fn'k-Q 



w - e„'k-Q - WqA - iO+ 

Nqx (T) + /„,k-Q 



(Bl) 



to - en'k-Q + c^qA - ^0^ 
In Ea lBll the dq integral is calculated using a numerical 



Montecarlo technique and taking explicitly into account 
the q — ?► divergence of |5^„/kp: 



dq- 



n 



RL 






dqcC 



(B2) 



In Eq |B2l the three-dimensional q integration compen- 
sates the q~^ divergence making the numerical evalua- 
tion of Eq|BT]feasible. Moreover, while | g„„/i5. | diverges 
as|Q|^-^,|Q|^| g^n'-k I is regular when Q — > 0. 
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